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ABSTRACT 

The evolution of present-day fossil galaxy groups is studied in the Millennium Simu- 
lation. Using the corresponding Millennium gas simulation and semi-analytic galaxy cata- 
logues, we select fossil groups at redshift zero according to the conventional observational 
criteria, and trace the haloes corresponding to these groups backwards in time, extracting the 
associated dark matter, gas and galaxy properties. The space density of the fossils from this 
study is remarkably close to the observed estimates and various possibilities for the remaining 
discrepancy are discussed. The fraction of X-ray bright systems which are fossils appears to 
be in reasonable agreement with observation, and the simulations predict that fossil systems 
will be found in significant numbers (3-4% of the population) even in quite rich clusters. We 
find that fossils assemble a higher fraction of their mass at high redshift, compared to non- 
fossil groups, with the ratio of the currently assembled halo mass to final mass, at any epoch, 
being about 10 to 20% higher for fossils. This supports the paradigm whereby fossils repre- 
sent undisturbed, early-forming systems in which large galaxies have merged to form a single 
dominant elliptical. 

Key words: cosmology: theory — galaxies: formation — galaxies: kinematics and dynamics 
— hydrodynamics — methods: numerical 
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1 INTRODUCTION 

Galaxy groups are believed to play a key role in the forma- 
tion and evolution of structure in the universe as, within a hi- 
erarchical framework, they span the regime between individ- 
ual galaxies and massive clusters. They are also more varied in 
their properties than galaxy clusters, as seen when various scal- 
ing re l ations are 
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compared with those of galaxy clusters jKaisei 
faD '199^; ' White et alj ri997|- lAUen & Fabiar 
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Xuec&Wu 2000; Helsdon & Ponmanl l200o[ 
Helsdon & Ponman 2003) . For instance, the relation between the 
luminosity and temperature of the X-ray emitting hot intergalactic 
medium (the L—T relation) has a larger scatter and a different slope 
for groups, when compared to similar properties of clusters. Vari- 
ous feedback mechanisms are often invoked to explain these dif- 
ferences. In addition, due to their lower velocity dispersion, groups 
are rapidly evolving systems, and galaxy mergers within groups can 
have a more significant effect on these relations than in clusters. 
In principle, the presence of cool cores and active galactic nuclei 
(AGN), as well as the star formation history, are all affected by ma- 
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jor interactions in the heart of a group or cluster. It would therefore 
be useful to find a class of groups or clusters with no major mergers 
in their recent history, to provide a baseline for the evolution of a 
passive system, with no major disruption. 

Fossil groups are good candidates for such a class of objects. 
They are distinguished by a large gap between the brightest galaxy 
and the fainter member s , wi th an under-abundance of L, galaxies. 
IZabludoff&MuIchaevI l Il998i) suggest that, for an X-ray detected 
group, the merging timescale for the most luminous group mem- 
bers {L ~ L*) is of order of a few tenths of an Hubble time, in 
agreement with the numerical simulations. A single giant elliptical 
;alaxy can fo rm as a result of multiple mergers within a few Gyr 
Barnej[T98^ . Thus, it is likely that one can find merged groups 
in the form of an isolated giant elliptical galaxy with an extended 
halo of hot gas, since the timescale for gas infall is longer than that 
on which galaxies merge (Ponman & Bertram 1993). In such sys- 
tems, the brighter galaxies, which have a relatively shorter merging 
timescale, are expected to me rge earli er leaving the fainter end of 
the luminosity function intact jPubins ki 1998; Miles et al.ll2004) . 

Following the discovery of a fossil gr oup having the abov e 
characteristics from ROSAT observations ^Ponman et al.l [l994h. 
more f ossil systems have bee n identified ( [Mulchaev & Zabludofl 
1 19991 ; IVikhlinin et 1 19991 ; Ijones et al.l 12000 ; IRomeii bOOoF 
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lMatsushitall200ll : |jones et alj|2003l). They are genera lly based on 
the definition of fossil groups from ljones et alj ( l2003h . i.e. groups 
with a minimum X-ray luminosity of Lx.boi ~ 0.25 x 10*'^/i~^erg 
s~^, as well as minimum magnitude difference of two between 
the first and second ranked galaxies, within half the pro- 
jected radius that encloses an overdensity of 200 times the 
mean density of the u n iverse (7?20o)- For an NFW profile 
( iNavarro. Frenk & White! Il996h . this is roughly equivalent 
to -Rsoo, the radius enclosing an overdensity of 500 times 
the mean (for NFW haloes of the appropriate concentration, 
Jisoo ~ 0.59 X R200). A few of these fossil groups have been the 
subje ct of detailed inve s tigatio ns ([Khosroshahi, Jones & Ponman 
2004 lYoshioka et id] |2004 [SunetalJ 120041: lUlmer et^ 
20051 : ICvpriano, Mendes de Oliyeira & Sodr j I2OO6I: 



Mendes de Oliveira. Cvpriano & Sodrd l200'i l Khosroshahi et al.l 
20061) . 



While most previous studies have focused on X-ray prop- 
erties of fossils, there is also emerging evidence that the galaxy 

S iroperties in fossils are different from those in non-fossils 
Khosroshahi. Ponman & Jones! [200^ . For instance the isophotal 
shapes of the central fossil galaxies appear to be non-boxy, sug- 
gesting that they may have formed in gas rich mergers. Vari- 
ous observational and theoretical studies ha ve suggested a signif- 
icant fractio n of galaxy g roups to be fossils 1 Vikhlinin et alJ 199S ; 
Jones etalj|200 3: D'Ong hia et ai]!2005! : lMilosavlievic et dxhmt : 



Sommer-Larsen 2006), though often the criteria used to define fos- 



sils in theoretical work are not easy to relate to observational stud- 
ies. 

Fossils may represent extreme examples of a continuum of 
group properties ~ they are consistently found to be outliers in 
the usual scali ng relations involving optical. X-ray and dynam- 
ical properties ^Khosroshahi. Ponman & Jones! !2007i) . While fos- 
sils fall on the L-T relation of non-fossil groups and clusters, 
they appear to be both hotter and more X-ray luminous than 
non-fossils of the same mass. Cooler fossil groups also show 
lower entropy than their non - fossil counterparts. According to 
!Khosroshahi. Ponman & Jones! j2007!) . the haloes of fossil groups 
appear to be more concentrated than those of non-fossil systems, 
for a given mass, which suggests that fossils have an early forma- 
tion epoch. As such, we have much to learn from them, and the in- 
vestigation of objects with similar properties in cosmological sim- 
ulations can provide important insights into the physical processes 
that underly the scaling relations. It can also reveal limitations in 
the numerical simulations, related to the treatment of physical ef- 
fects like pre-heating, feedback and merging, which are difficult to 
model. It is thus important to study the formation and evolution of 
these systems in the cosmological N-Body simulations which have 
become essential tools for studying formation of large scale struc- 
ture in the Universe. 

^In this paper we use the Millennium simul ation 1 Springel et al.! 

!2005l) together with the semi-analytic models jCroton et alJ!200^ 
of galaxy form ation within dark m atter haloes and the Millennium 
gas simulation jPearce etal ]|2Q03), to identify fossil groups, study 
their properties in the simulations and make a comparison to the 
observations. We begin with a brief discussion in §2 of the Millen- 
nium Simulation, and the implemented semi-analytic galaxy cata- 
logues and gas simulations. In §3 we discuss our method of identi- 
fying optical and X-ray fossil groups from these catalogues. In §4, 
we discuss the various properties of these fossil groups, their abun- 
dance in the local Universe and the evolution of simulated X-ray 
fossils with time. Finally, in §5, we summarize the implications of 
our results in terms of the evolution of fossil groups in the context 



of multiwavelength observations. Throughout the paper we adopt 
Ho = 100 h km s~^ Mpc~^ for the Hubble constant. 



2 DESCRIPTION OF THE SIMULATIONS 

In order to extract fossil groups in the Millennium simulation, us- 
ing observational selection criteria, we require a simulation suite 
that includes the baryonic physics of hot gas and galaxies, as well 
as a high resolution dark matter framework and a sufficient spatial 
volume to limit the effects of cosmic variance. Fo r this study we use 
the dark matter Millennium Simulation ( Springel et al. I2OO5I) , a 10- 
billion particle model of a comoving volume of side 500/i~^ Mpc, 
on top of which a publicly available semi-analytic galaxy model 
teroton et all!200 ^ has been constructed. For the hot gas we have 
repeated the Millennium simulation with a lower resolution simu- 
lation including gas physics utilising the same volume, phases and 
amplitudes as the original dark-matter-only model. This run accu- 
rately r eproduces the struc tural framework of the Millennium Sim- 
ulation jPearce et al.l2007l) . Below we summarize the main charac- 
teristics of the above simulations. 



2.1 The Millennium Simulation 

The Millennium Simulation is based on a Cold Dark Matter cosmo- 
logical model of structure formation, with a Dark Energy field A. 
The basic assumptions are those of an inflationary universe, dom- 
inated by dark matter particles, leading to a bottom-up hierarchy 
of structure formation, via collapsing and merging of small dense 
haloes at high redshifts, into the large virialised systems such as 
groups and clusters that contain the galaxies that we observe today. 
The simulation was performed using the publicly available paral- 
lel TreePM code Gadget2 ( Springel et al. 2001), achieving a 3D 
dynamic range of 10''' by evolving 2160"^ particles of individual 
mass 8.6 x lO^h^^ Mq, within a co-moving periodic box of side 
5QQh~^ Mpc, and employing a gravitational softening of 5h~^ kpc, 
from redshift z = 127 to the present day. The cosmological pa- 
rameters for the Millennium Simulation were: f^A = 0.75, Q.m = 
0.25,^6 = 0.045, /i = 0.73,71 = 1, and erg = 0.9, where the 
Hubble constant is characterised as 100 /i kms~^Mpc~^. These 
cosmological paramet ers are consistent wi th recent combined anal- 
ysis from WMAP data (Spergel et alj2003l) and the 2dF galaxy red- 
shift survey i Colless et al. .2001„) , although the value for ag, is a little 
higher than would perhaps have been desirable in retrospect. 

The derived dark matter halo catalogues include haloes down 
to a resolution limit of 20 particles, which yields a minimum halo 
mass of 1 .72 x 10^" ^ Mq . Haloes in the simulation are found us- 
ing a friends-of-friends (FOF) group finder, tuned to extract haloes 
with overdensities of at least 200 relative to the critical density. 
Within a FOF halo, substructures or subhaloes are identified using 
the SUBFIND algorithm developed by Springel et al. (2001), and 
the treatment of the orbital decay of satellites is described in the 
next section. 

During the Millennium Simulation, 64 time-slices of the loca- 
tions and velocities of all the particles were stored, spread approx- 
imately logarithmically in time between z = 127 and 2 = 0. From 
these time-slices, merger trees are built by combining the tables of 
all haloes found at any given output time, a process which enables 
us to trace the growth of haloes and their subhaloes through time 
within the simulation. 
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2.2 The Semi Analytic model 

Using the dark mat t er hal oes of the ISpringel et all bOOSi) simu- 
lation, ICroton et alj ( l2006h have simulated the growth of galax- 
ies, and their central supermassive black holes, by self-consistently 
implementing semi-analytic models of galaxies on the outputs of 
the Millennium Simulation. The semi-analytic catalogue contains 
9 million galaxies at 2 = down to a limiting absolute magnitude 
of A/_R — 5 log h = —16.6, observed in B, V, R, I and K filters. 
The models focus on the growth of black holes and AGNs as feed- 
back sources. The inclusion of AGN feedback in the semi-analytic 
model (which allows the cooling flow to be suppressed in massive 
haloes that undergo quasi-static cooling) and its good agreement 
with the observed galaxy luminosity function, colour distribution 
and the clustering properties of galaxies, make this catalogue well- 
matched and suitable for our study of fossil systems. 

In the semi-analytic formulation, galaxies initially form within 
small dark matter haloes. As the simulation evolves, such a halo 
may fall into a larger halo. The semi-analytic galaxy within this 
halo then becomes a satellite galaxy within the main halo and fol- 
lows the track of its original dark matter halo (now a subhalo) until 
the mass of the subhalo drops below 1.72x lO^^/i^^ M0, which 
corresponds to a 20-particle limit in the Millennium Simulation. At 
this point the galaxy is assumed to spiral into the centre, on some 
fraction of the dynamical friction ti mescale, wh ere it merges with 
the central galaxy of the larger halo jCroton et alj2006i) . 



2.3 The Millennium Gas Simulation 

The Millennium Gas Simulations are a suite of hydrodynamical 
models, utilising the same volume, and values of initial perturbation 
amplitudes and phases as the parent dark-matter-only Millennium 
Simulation. Each of the three models completed to date contains 
additional baryonic physics: The first does not follow the effects 
of radiative cooling and so overpredicts the luminosities of group- 
scale objects significantly. The second includes a simple preheat- 
ing scheme that is tuned to match the observed X-ray properties of 
clusters at the present day and the third includes a simple feedback 
model that matches the observed properties of clusters today, as 
well as having some chance of following the time evolution. We 
have used the second of these models in this work, as we only 
utilise the hydrodynamical properties of the groups at 2: = 0, where 
the observational and simulation results are well matched. 

Each of the Millennium Gas Simulations consists of 5 x 10* 
particles of each species, resulting in a dark matter mass of 1.422 x 
10^°h~^ Mq per particle and a gas mass of 3.12 x lO'^h"^ M© 
per particle. The Millennium Simulation has roughly 20 times bet- 
ter mass resolution than this and so some perturbation of the dark 
matter halo locations is to be expected. In practice the position and 
mass of dark matter haloes above 10^^ Mq are recovered to 
within 50 /i~^kpc between the two volumes, allowing straightfor- 
ward halo-halo matching in the large majority of cases. 

The Millennium gas simulations used exactly the same cos- 
mological parameters as those stated above. With the inclusion of a 
gaseous component, additional care needs to be taken in choos- 
ing the gravitational softening len gth in order to avoid spurious 
heating jSteirmietz & Whitall997l) . We use a comoving value of 
25(1 -1- z)h~^ kpc, roug hly 4% of the mean interparticle separa- 
tion teorgani et al.ll200q ) until 2 = 3, above which a maximum 
comoving value of 100/i~^ kpc pertains. We have adopted a differ- 
ent output strategy for the Millennium Gas Simulations, preferring 
to output uniformly in time with an interval roughly corresponding 
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Figure 1. The bolometric X-ray luminosity versus dark matter halo tem- 
perature (left panel) and mass within -R200 (right panel) for all haloes 
in the Millennium gas simulation. The vertical dashed lines correspond 
to the X-ray luminosity threshold Lx.bol = 0.25 X lO*'^ ^"^erg s"^ 
adopted in this paper for defining fossil groups. The cutoff M(iJ2oo) ^ 
IQi^ fc-l Mq is adopted in fi \3l\ 



to the dynamical time of objects of interest. This strategy results in 
160 rather than 64 outputs and places particular emphasis on the 
late stages of the simulation. 



3 SAMPLE SELECTION 
3.1 Definition of fossils 

Fossil groups are selected according to a combination of X-ray 
and optical criteria, based on the observational definition given by 
I Jones et ^ ( |2003 \ which is widely followed in the literature. Their 
X-ray luminosity must satisfy Lx.boi ^ 0.25 X lO'^'^h-^ erg s , 
and the difference between the i?-band magnitudes of the first and 
second ranked galaxies, within half the projected radius enclosing 
200 times the mean density (7?20o), must be Ami2 ^ 2 magni- 
tudes. A limit of 0.5i?200 is used because i* galaxies within this 
radius should spiral into the centre of the grou p due to orbital de - 
cay by dynamical friction within a Hubble time jjones et alj2003h . 
The limit on the bolometric X-ray luminosity Lx.boi helps to ex- 
clude poor groups and individual galaxies with a few small satel- 
lites. Such groups are often not in dynamical equilibrium, and in 
addition there might be a gap in their galaxy luminosity function 
simply as a result of the small numbers of galaxies involved. We 
address this issue in 13.31 



3.2 Selection of Fossil galaxies 

For this investigation, we first selected dark matter haloes 
from the Millennium Gas Simulation with masses M{R2oo) ^ 
10^^ Mq. As Fig. [T] demonstrates, all haloes for which 
Lx,boi> 0.25x10^2 erg s ^ are expected to be included in this 
sample and thus our fossil sample will be complete. A bolometric 
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X-ray luminosity of 0.25 x 10*^ erg s^^ corresponds to a tem- 
perature of T~0.5 keV. This constraint on tlie halo masses results 
in an initial sample of 51538 dark matter haloes, within which we 
search for fossil groups. 

In order to extract information about dark matter and galaxy 
properties, for each of these 51538 haloes found in the Millen- 
nium Gas Simulation, the counterpart in the Millennium Simula- 
tion needs to be found. This is a straightforward procedure because 
the simulated volume and the amplitudes and phases of the ini- 
tial power spectrum were matched. However, the dark-matter-only 
Millennium Simulation has 20 times the mass resolution of the cor- 
responding Millennium Gas Simulation and so contains some ad- 
ditional small scale power. This leads to small offsets (typically 
around twice the gravitational softening length) in the final co- 
ordinates of equivalent haloes, and even smaller mass differences 
(typically less than 5%). Of the 51538 haloes, 48774 (~95%) have 
corresponding haloes identified in the Millennium Simulation. For 
each of these matched haloes we extracted the coordinates and 
BVRIK magnitudes for each gala xy contained within R 200, from 
the publically available catalogue of lcroton et all ilOOdtl . 

The simulated properties of the galaxies occupying each dark 
matter halo were then used to calculate Ami2, i.e. the differ- 
ence in i?-band magnitude of the first and second ranked galaxies 
within 0.5i?200 of the centre of the halo. Out of the 48774 matched 
haloes, 6502 are found to be optical fossil groups, i.e. haloes with 
Ami2 2 mag, among which 1300 are X-ray fossil groups, i.e. 
optical fossil groups with Lx,boi ^ 0.25 x 10*^ erg s~^. As 
can be seen from Fig.[2l X-ray fossil groups do not form a separate 
population but are rather extreme examples of a smooth distribu- 
tion. It is also clear that the spread in Ami2 increases dramatically 
with decreasing in the X-ray luminosity and hence the enclosed 
mass, in agre ement with the condi t ional l uminosity function (CLP) 
formalism of Ivan den Bosch et al.l j2007l) . 

Control groups are necessary for our study in order to allow 
us to compare the properties of X-ray fossil and non-fossil groups. 
We define two control group samples, based on the magnitude dif- 
ference of the two brightest members of the group (within 0.5-R200 
of the centre of the dark matter halo): (i) 0.8 ^ Ami2 5S 1.0, and 
(ii) 0.1 ^ Ami2 ^ 0.3. While members of the former are examples 
of intermediate groups, the latter could be regarded as a class of ex- 
treme non-fossil groups in that they contain at least two galaxies of 
very similar magnitude. The bolometric X-ray luminosity limit for 
each of the control samples is the same as that of the X-ray fossil 
group sample (see Fig.|2l(. 

3.3 The likelihood of finding groups with Ami2 ^ 2 at 
random 

One of the central criteria used to define fossils is the absence of 
galaxies within a range of two magnitudes of the brightest galaxy 
(Ami2 ^ 2). However, for groups with only a small number of 
members, there is a significant probability of obtaining such a lu- 
minosity gap as a natural consequence of the high-end tail of the 
galaxy luminosity distribution. To quanti fy the likelihood o f ob- 
taining a value of Ami2 ^ 2 by chance, Ijones et ^ ( |2003|) per- 
formed 10* Monte Carlo simulations for groups and clusters with 
abso lute magnitudes selected at random from a Schechter function 
jSch echter 1976). Using the paramet ers of the composite luminos - 
ity function of MKW/AWM clusters ^Yamagata & Maeharalll986h , 
they found that for the systems of ~40 galaxies, 0.4±0.06% of the 
generated luminosity functions had Ami2 2. 

We performed a similar analysis for groups spanning a range 




^m^, (within O.SR,^^^) % of lialoes {Mr\^>2) 

Figure 2. X-ray luminosity versus the i?-band luminosity gap Ami 2 
within 0.5i?.2oo ench of the dark matter haloes with gas properties. The 
horizontal dashed-line intersects the vertical axis at Lxfiol = 0.25 X 
lO'^^h^^erg s^^. The top-right part of the graph shows the region for 
which Ami2 ^ 2 mag and Lx,bol ^ 0.25 X 10*^/i~^erg s~^, the 
optical and X-ray criteria that jointly define fossil groups. The two shaded 
regions show the location of our control samples (see >i l3.2l for details). The 
histogram shows the fraction of optical fossils in each bin of Lx^bol- 



in richness, and using parameters appropriate to our data from 
the Millennium simulation. For twenty classes of groups, contain- 
ing 10, 15, 20, 25,... galaxies respectively, we randomly generated 
galaxies according to a Schechter luminosity function. (None of our 
X-ray fossil groups from the Millennium simulation contain fewer 
than 10 galaxies.) The characteristic magnitude M* ~ —22.1 and 
faint-end exponent a ~ —1.19, were adopted from a fit of the 
Schechter function to the i?-band luminosity function of the semi- 
analytic catalogue. A magnitude cut off of —17.4 was then applied 
to the magnitude of generated galaxies as this is the i?-band mag- 
nitude completeness limit of the semi-analytic catalogue. lO'^ sim- 
ulations were carried out for each richness class of group. 

Fig. [3] compares the percentage of optical and X-ray fossil 
groups from the Millennium simulation as a function of number 
of galaxies within 0.5i?200 for each dark matter halo, with those 
populated using a Schechter luminosity function as detailed above. 
The lower panel shows the result when the expected number of ran- 
domly generated groups with Ami2 ^ 2, is subtracted from the op- 
tical fossils. For poor systems, the incidence of 'statistical fossils' 
is significant. It can be seen that approximately one third of the 
fossil systems with fewer than 25 galaxies seen in the Millennium 
data can be attributed to statistical chance (as opposed to the result 
of physical processes generating a non-statistical luminosity gap). 
However, even after these random fossils are removed, the fraction 
of optical fossils increases as the number of galaxies within dark 
matter haloes decreases. In contrast, for X-ray fossils, many of the 
poor haloes which qualify as 'statistical fossils' fail to pass the X- 
ray luminosity threshold criterion, so the chance fraction is never 
much larger than 20%. We return to this issue in Sec. 14.1.21 For 
groups with more than 30 galaxies, the fraction of fossils meeting 
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Figure 3. The histogram shows the incidence rate of Ami2 ^ 2 occurring 
by chance from a random population of galaxies selected from a Schechter 
luminosity function in comparison to fraction of optical fossils (open cir- 
cles) and X-ray fossils (solid triangles) in the Millennium simulation, as a 
function of number of galaxies per halo within 0.5-R200 ■ The lower panel 
plot (filled circles) shows the result of subtracting random groups from the 
optical fossil groups. 



the Ami2 ^ 2 criterion by chance drops below 1.0%, and soon 
becomes negligible. 



Figure 4. The R-band luminosity gap distribution for haloes from the Mil- 
lennium semi-analytic model within the mass range M = 0.5 — 10 X 
10^* Mq, evaluated relative to the first and second most luminous 
galaxies and the first an d the third most luminous galaxies, superposed on 
the theoretical model of iMilosavlievic et alj j200 6) and the SDSS data for 
the same mass range but different searching radius (-R200 for the model 
and projected radius of 500h~^kpc for the SDSS data). The Millennium 
data are plotted within -R200 (closed circles) as well as the projected radius 
of 500/i~^kpc (triangles), (a), (b): The lumino sity gap statistic predictions 
of the theoretical model of lMilosavlievic et al.l I2OO6.) with In A = 1 (thin 
green line) and InA = 2 (thick blue line), (c), (d): The r-band luminosity 
gap dist ribution from 730 c lusters (red histogram) in the SDSS C4 Cata- 
logue of lMiller et all )2005l) . 



4 RESULTS 

4.1 The Luminosity Gap Statistic 

If galaxies in groups and clusters merge with the central galaxy 
over a finite time to produce a progressively greater massive cen- 
tral galaxy, then one way to quantify the dynamical age of a galaxy 
system is to measure the luminosity gap between the two most lu- 
minous galaxies remaining in the group, provided that there is no 
infall of bright galaxies from other nearby systems. A relaxed X- 
ray morphology in observed fossils is an indication of the absence 
of such a process. Within such a merging scenario for the forma- 
tion of the brightest group galaxy, the luminosity gap distribution 
depends on the halo merger ra tes. This has pre viously only been 
examined in analytical studies l lBond et al.|[T99ll) , or in large-scale 
cosmological simulations without a hot baryonic component. Thus, 
while an important observational criterion for defining fossils relies 
on the X-ray properties of the group, this has not previously been 
implemented within the models. 



4.1.1 The R-band magnitude gap distribution 

Here we compare the luminosity gap distribution between the lu- 
minous galaxies found at the centre of dark matter haloes extracted 



from the Millennium simulation with expectations from the ana- 
lytical model of Milosavljevic et al,. (2006) , and with observational 
prop erties from SDSS clusters. 

IMilosavlievic et al] 1 I2OO6I) compared the distribution of the 
predicted luminosity gaps from their analytical model within -R200, 
as a function of halo mas s, with the observ ed luminosity gaps in 
the SDSS (DR4) clusters tMiller et alj|2005l) . ranging in mass and 
redshift from M = 0.5-10 x lO" Mq and 2=0.02 to 0. 17, re- 
spectively, within a projected physical radius of 500 /i~^kpc. Halo 
merger rates in their model have been analytically estim ated ac- 
cording to the excursion-set theory of iBond et al.l ( Il99 ll) . which 
is also known as the extended Press-Sche chter formalism. Assum- 
ing a halo density profile of the form of iNavarro. Frenk & White! 
( Il996r) . a subhalo of mass m merges into a primary halo of mass M 
(m ^ ^M) and makes a composite halo. As the centre of the sub- 
halo crosses the virial radius of the new composite halo, a merger 
happens. Then the subhalo spirals toward the centre of the compos- 
ite halo in a near circular orbit, experiencing dynamical friction. 

For comparison with the above study, we evaluate the _R-band 
luminosity difference between the first and second most luminous 
(Ami2) and the first and the third most luminous (Amis) galax- 
ies in our Millennium data within -R200 and within 500 h^^kpc. 
In Fig. |4] we plot the i?-band luminosity gap distribution of the 
Millennium simulation for the same mass range as the models, to- 
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geth er with the lumino sity gap distribution of 730 SDSS C4 clus- 
ters jMilleretal .1120051). Figs. [4^1 and c ompare the predicted 
gap statistics from' MilosavlieviJ et al.l l l2006h for two values of the 
Coulomb logarithm, In A = 1 and In A = 2, within -R200 • Since 
the parameter In A is proportional to the force of dynamical fric- 
tion between the centres of subhalo and primary halo during the 
process of merging, a higher value of In A corresponds to a faster 
effective halo merger rate. In numerical simulations, In A is ap- 
proximated by bmax/bmin, whcrc bmax and bmin are the maxi- 
mum and minimum imp act parameters respectively, and In A is ex- 
pected to be ~ 1 - 4 jVelazguez & White|[T99a: iFellhauer et~al] 
[2000 ; D'Ongh ia et alj20q5l). However, in the semi-analytic galaxy 
catalogues ( Crotonet al.||2006l) , based on the Millennium simu- 
lation, used in this work, the above relation is approximated by 
InA = ln(l + M200 / msax) , where nisat is the halo mass of the 
satellite galaxy. 

Within the mass range of the SDSS data there are 8842 haloes 
in the Millennium simulation catalogue. Accordingly, in Fig.lD our 
data have been normalised t o be comparable with the SDSS data 
and the theoretical model of iMilosavlievic et al.l ( l2006h . However, 
the simulation data is, unlike the observations, complete and un- 
contaminated by spurious groups or foreground and background 
galaxies. All these effects are likely to be heavily dependent on the 
number of galaxies residing in the halo. As such, the comparison 
with the SDSS data shown in Fig.|4]should be treated with caution. 

Given this caveat, our analysis based on the Millennium sim- 
ulation catalogues agrees remarkably well with the models of 
IMilosavlievic et alj (2006) based on the SDSS survey for the lu- 
minosity gap distribution of the two brightest galaxies in each of 
the dark matter haloes, particularly for InA = 2 (Fig. |4ji). How- 
ever, for the _R-band luminosity gap between the brightest and third 
brightest galaxies in each system (Fig.|4j5), the simulations signif- 
icantly depart from the model. When comparing with the SDSS 
data (Figs.|4j; and |4jl), the simulations overpredict the frequency 
of groups. The simulations and the SDSS data have similar shaped 
distributions for the luminosity gap Amis, but with a shift of ~0.5 
mag toward higher Amis in the simulated haloes. 

We emphasize that the Millennium predictions for the lumi- 
nosity gap statistic are sensitive to the assumed mass range and 
search radius of dark matter haloes within which brightest halo 
members are identified. SDSS cluster masses have been estimated 
from total r-band luminosities, so any inaccuracies in this proce- 
dure would affect the comparison with the Millennium data. 

Observationally there is an excess population of groups with 
a small luminosity gap between the first and second ranked galax- 
ies, above what is predicted by the theoretical models or the sim- 
ulations. This excess population is likely to result from contam- 
ination of observed group samples by local structure alignments, 
and renormalising to a sample without these groups scales down 
the "Millennium" distribution in Fig. |4j;, bringing the simulation 
results and the observational measurements into better agreement. 
Results are similar in the if-band. 



4.1.2 The abundance of fossil groups 

The probability of finding fossil systems is expected to in- 
crease with decreasing halo mass, as shown in previous stud- 
ies based on theoretical mod els or hydrodynamical simulations 
(D'Onghia et al. 2005; Milosa vlievic et alj|20o3 : ISommer-LarsenI 
boofe .van den Bosch et al..,2007D . Unfortunately, it is difficult to 
compare the results from different studies (both theoretical and ob- 
servational), since they have used a range of search radii (from 



J?i8o to i?337 - see Table. 1) within which the Ami2 ^ 2 mag 
criterion is imposed. Clearly, the larger the search radius, the more 
demanding is the requirement on the galaxy contents of the system, 
and the smaller the fraction of groups which will qualify as fossils. 

In Fig.|5] the rates of incidence, Pf{M), of optical fossils and 
X-ray fossils (using our preferred search radius of 0.5i?2oo, fol- 
lowing Oones et al.l ( l2003h ) are plotted, as a function of the mass M 
of the halo, togethe r with the predicted values from the models of 
IMilosavlievic et al] ( I2OO6I) for two values of A. The shape of our 
curve for optical fossils is quite similar to the theoretical models 
(which included no X-ray luminosity criterion), but the latter ac- 
tually employed a search radius of J?200- To see the effect of this, 
we also show our Millennium results for this larger search radius. 
The fraction of fossil systems falls by approximately a factor of 
2, when this more demanding requi rement is imposed, and so lies 
significantly below that predicted bv lMilosavIievic et Sllzooe"). 

On scales of M ~ 10^^ - 10^*/i"^Mq, ~5%-18% of groups 
are optical fossils. This probability falls to ~3%-5% for more 
massive (M ^ IO^'^/i^^Mq) fossil systems. For halo masses 
> 5 X 10^"^ h^^MQ all optical fossils in the simulation are also 
X-ray fossils. However, at the lowest halo masses the fraction of 
X-ray fossils drops steeply, since many low mass haloes do not sat- 
isfy the Lx threshold criterion. 

In Table. 1, we summarize the incidence rates of fossil systems 
from present study as well as those found in the literature. Compar- 
ison between these different estimates is difficult, since both the 
search radius and the halo mass range varies considerably from 
study to study. However , a direct comparis on with the only obser- 
vational estimate (from ljones et al.l l l2003h ) is possible, since we 
have used the same definitions of fossil groups as these authors. 
Based on a comparison w ith the integrated local X -ray luminosity 
function of .Ebeling et alj ( l200ll) . |jones et all ( l2003h estimated that 
X-ray fossil systems constitute 8-20% of all systems of the same X- 
ray luminosity (Lx.boi ^ 0.25 x 10*'^/i~^erg s~^). The right panel 
histogram of Fig.|2|represents the fraction of optical fossil systems 
in each bin of Lx,boi- Integrating this over all X-ray luminosities 
above the threshold value for fossils, we find that ~ 7.2 ± 0.2% 
of haloes with Lx,boi ^ 0.25 x 10^^/i^^erg s^^ are X-ray fos- 
sils, which is reasonably consistent with the lower limit of ~ 8%, 
derived by Jones et al. (2003.). 

In comparison, detai led hydrodynam i cal s imulations by 
iD'Onghia etal] | |2005|) and ISomnier-Larsen! 1 I2OO6I) of 12 galaxy 
groups, predict a larger fraction of 33%±I6% for fossil systems 
of mass W^'^Ii'^Mq or larger. This may be because it is easy 
to overestimate th e local viscosity in hydrodynamic simulations 
jTittlev et aLll200lh . a process that would lead to central overmerg- 
ing in the models. 



4.2 The Space Density of X-ray Fossil Groups 

So far, the integrated space density of X-ray fossil groups has 
been studied for small samples, each of th ree to five X-ray fos- 
sil systems, at different limiting luminosities jVikhlinin etaljl999l; 
lRomeill200d Ijones et al.ll2003[) . Here, we estimate the space den- 
sity by systematically counting the fossil groups in the whole 500 
Mpc survey volume of the Millennium Simulation at z = 0. 
For comparison with previous studies, we select and count X-ray 
fossil groups for three limiting X-ray luminosities ranging from 
0.25-5x10''^ fe-^erg s-\ The space densities calculated at differ- 
ent limiting luminosities as well as those from previous studies are 
given in Table 2. The value from Romer et al is a very rough esti- 
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Table 1. The incidence rates of fossil systems. 



Mass range (/i ^Mq) 


Lx (10*2 ergs-i) 


Fossil type^ 


Search radius 


Fossil fraction (%) 


Reference_^ 






O 


R200 


~ 5 - 40 


M06 


~ 1013 -10" 




O 


RlSO 


~ 3.6 ±0.1 


vdBOT^ 


> 10" 




o 


R200 


~ 1 - 3 


M06 ~ 


~ 10" -1015 




o 


R18O 


~ 6.5 ±0.1 


vdB07 


-10" 




o 


^?337 


~ 33 ± 16 


SL06, DOOsf, 


~10i3-10i5 




o 


1 fc-lMpc 


~ 8 - 10 


S07^ ~ 




^ 0.25 


X 


0.5iJ2O0 


~ 8 - 20 


J03~ 




^ 0.25 


X 


0.5iJ200 


~ 7.2 ±0.2 


Present study-^ 


~ 1013 -1015 




o 


0.5i?200 


~ 13.3 ±0.2 


Present study 



" O: Optical fossils , X: X-ray fossils. 

^ S07 : ISales et al.l fcOOTh: M06: iMilosavlievic et alj j200d): vdBOV: Ivan den Bosch et alj <2007l) : SL06: ISommer-LarsenI 
j20O6l) : DO05: lD'Onghia et"ai] i2005l) : J03: |jones et alj j2003l) . 

Based on the conditional luminosity function (CLF) formalism of lvan den Bosch et alj i2007l) . 

From hydrodynamical simulations of 12 galaxy groups. 
^ Based on the Millennium simulation. The first brightest galaxies of fossils in their sample are always brighter than Mn = 
-20.5. 

f Histogram on the right panel of Fig.|2] gives the fraction of X-ray and optical fossils in each bin of Lx ■ 



Table 2. Space densities of fossil galaxy groups. 



— Ln A=1; Rj, 

Ln A =2; R^, 

• Optical fossils; R^,, 
O Optical fossils; 0.5R^^ 
X-ray fossils; 0.5R 




N, 



Density_ 



Reference_ Present study'' 



13.5 14.0 14.5 

Log (M/h"'M . ) 



Figure 5. The probability, Pf (M), that a dark matter halo of mass AI con- 
tains an X-ray fossil group (closed triangles), optical fossil group within 
0.5i?200 (open circles), or optical fossil group within i?200 (closed circles) 
from the Millennium simulation. The fossil incidence rate from the analyt- 
ical study of Milosavljevic et al. 1 2006) for two values of LnA = 1 (green 
line) and LnA = 2 (blue line) is also plotted. The vertical dotted-Hne cor- 
responds to halo mass ~ 3.34 X IOI3 fc-l M0 (see Sec.l431. 



mate, since no redshifts for galaxies surrounding the central object 
were available in this study. 

Our values show that for X-ray luminosities exceeding 2.5- 
SxlO'*'^ h~^eTg s~i, the space density of fossils in the Millen- 
nium simulation agrees within the errors with those estimated by 



> 0.25 


5 


320l?« 


J03 


104 ±3 


> 2.5 


3 




J03 


22.4 ± 1.3 


> 2.5 


4 




V99 


22.4 ± 1.3 


> 2.5 


3 


- 160 


ROO 


22.4 ± 1.3 


> 5.0 


4 


1 q 2+2'i-** 


V99 


12.8 ± 1.0 



In units of lO^^/i'^ erg s"! 
^ Number of fossils 

In units of IQ-'^h^ Mpc'^ 

V99 lVikhlinm et aljh999l) . ROO lRomeJ ilOOdi . J03 [jones et aljjiooj) 



IVikhlininet"aLl ( ll999l) and |jones~et 'Zl( l2003l) . At the lowest X-ray 
luminosities, the density from the Millennium simulations appear 
to be lower than observed, though the observational values given 
in Table 2 have large uncertainties due to the small number of X- 
ra y fossil groups and the effects of cosmi c var iance. Recent studies 
of lKhosroshahi. Ponman & Jones! ( l2007h a ndljeltema et al. (200^ 
show that one of the fossils in the sample of l Jones et al.n2003) does 
not satisfy the fossil criterion of Am 12 ^ 2, which reduces the 
observational space density. Certainly the number of X-ray fossils 
found is heavily dependent on the X-ray luminosity threshold cho- 
sen and may be influenced by the scatter in X-ray group properties 
near this lower limit. 



4.3 Evolution of Fossil Groups 

Strong interactions and mergers between galaxies occur more effi- 
ciently in the low velocity dispersion environment of galaxy groups 
(Miles et al. 2004). Therefore in old, relatively isolated groups, 
most massive galaxies have sufficient time to merge via dynamical 
friction. If X-ray fossil groups are indeed systems that formed at 
an earlier epoch, we should be able to verify this from the merger 
histories of present-day fossils in the Millennium simulation: an 



8 Dariush et al. 



1.0 



0.66 



Redshift 

0.42 0.24 



0.11 



Redshift 

1.0 0.66 0.42 0.24 0.11 0.0 



T 



1.0-3 



0.9 



5 0.8 



0.7 - 



0.6 



0.5 - 



0.4 



it 

It ' 



1 ' 1 ' r 

i-^~2.8Gyr 



I 



• X-ray fossil groups 
A Control groups 1 

* Control groups 2 



,1' 



High-mass haloes 



H 1 h 



Low-mass haloes 



1.0 



0.8 



0.6 1 
3 

03 
cn 
m 

0.41! 

i.oi. 



0) 



0.5 0.6 0.7 0.8 0.9 

Expansion factor 



1.0 



0.5 0.6 0.7 0.8 0.9 1.0 

Expansion factor 



0.8 



0.6 



0.4 



Figure 6. Tracing back the mass build-up of the dark matter haloes as a function of expansion factor and redshift for both the X-ray fossils and the control 
groups 1 and 2. (a) For all haloes, (b) High mass haloes (m^^^i^ > 3.34 X 10^'' Mq). Both plots indicate the earlier formation of X-ray fossil groups in 
comparison to control groups, (c) Low mass haloes (mjjaio ^ 3.34 X lO'^^ Mq). Here the difference in evolution between the X-ray fossil and control 
groups is not as pronounced as in those seen in high mass X-ray fossils. All the masses are normalised to the mass at z=0. 



exercise that is not directly possible to perform with observational 
surveys. In Fig.|6|we trace the mass evolution of present-day X-ray 
fossil systems backwards from z = io z = 0.82 when the scale 
factor, a, of the Universe was 0.55 times its current size. 

At any given redshift in Fig. [6] the average ratio of the mass 
of a halo to its final mass (at z = 0) is calculated for all eligible 
haloes. The error is represented by the standard error on the mean, 
i.e. a/\/iV, where a is the standard deviation of the original distri- 
bution and A'^ is the sample size. The same was done for both sets 
of control groups. The original sample of fossils was divided into 
two subsamples, of low mass and high mass groups (see Fig.|6j5 and 
Fig.|6j;), such that both subsamples have equal numbers of groups. 
The boundary between the two subsamples corresponds to the me- 
dian present-day mass ~ 3.34 X 10" /i-^Mq. 

Fig. [6^ shows that at a scale factor of 0.8 {z ~ 0.24), the fossil 
groups have already attained ~90% of their final mass while, at the 
same redshift, the fraction of assembled mass of the extreme non- 
fossil groups is about ~77% of their final mass. The intermediate 
control group gives intermediate values. The fossil groups have al- 



most all their mass in place by a redshift of z ~ 0.1, and show no 
evidence of recent major mergers, while the non-fossils seem to be 
assembling mass even at the present day. These results suggest an 
early formation and consequent higher mass concentration in fossil 
groups, in comparison to normal groups, particularly for the more 
massive fossils. 

As Figs.|6j3 and [6}; show, the difference in mass assembly is 
larger in more massive haloes than haloes with lower mass. The 
decreased distinction in the assembly history for our lower mass 
fossil systems probably results from the fact there is a large frac- 
tion of "statistical fossils" in this category: groups which achieve 
Ami2 ^ 2 due to random chance, because of the small number of 
members. As can be seen in Fig.jH ~50% of optical fossil groups 
with masses less than ~ 3.3 x 10^^ M0 are expected to fall 
into this "statistical fossil" category. 



V ari ous obs erv ational 
1 1994 I Jones et alj l2003l : 

2004,; Ul mer et al.l [20051 

iKhosroshahi. Ponman & Jones! 



properties (jPonman et al 



IKhosroshahi . Jones & Ponmar 
[Khosrosha hi et alj l2006l: 
2OO7I) have suggested an early 
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formation epoch for fossils. iD'Onghia et alj ( |2005|) and 
ISommer-LarsenI ^20061) used a set of twelve high-resolution 
numerical simulations in the ACDM cosmology to study the for- 
mation of fossil groups, and found a correlation for the magnitude 
gap between the brightest and second-brightest galaxies and the 
halo formation epoch, with fossils accreting half of their final dark 
matter mass at z ^ 1. Such an early assembly of fossil haloes 
leaves enough time for L, galaxies to merge into the central one 
by dynamical friction, resulting in the observed magnitude gap at 

2 = 0. 

Fig.|7]shows the history of mass assembly of a typical example 
of a massive fossil group (right panel) and a control group (left 
panel) from the Millennium Gas Simulation from redshift z = 1.0 
to 2 = 0. The dimension of each image is 10 x 10 Mpc, centred 
around the central halo. It can be seen that al z — 0.3, the X-ray 
fossil group has already largely been assembled, while the control 
group has considerable substructure even at a later epoch. 



5 DISCUSSION 

We studied the history of the mass assembly of fossil groups, se- 
lected using the usual observational criteria at 2 = 0, from a red- 
shift 2 = 0.8 to the present day, within the Millennium simula- 
tion. A sample of X-ray fossil groups was defined from the Mil- 
lennium simulations and associated gas and galaxy catalogues, ac- 
cording to the usual criteria: (a) the difference between the i?-band 
magnitudes of the first and second ranked galaxies, within half the 
projected radius enclosing 200 times the mean density of material 
(^^200), is Ami2 2 magnitudes, and (b) The bolometric X-ray lu- 
minosity of the group is Lx.boi J? 0.25 x lO^^/i'^erg s~^. While 
optical fossil groups fulfill just the first condition. X-ray fossils sat- 
isfy both criteria. Our main results are as follows: 

• The space density of X-ray fossil groups is in close agreement 
with the observed space density of fossils with Lx > 2.5 x 10^^ 
h~^erg s~^. Although for low luminosity fossils we find roughly 
1/3 of the observed fossil space density, there are several poten- 
tial factors that could lead to this difference. As well as significant 
uncertainties in the observational studies, the X-ray properties of 
haloes in the real Universe show far gr eater scatter than tho se seen 
in the preheating simulation used here i Hardev et aklbOOTh . Given 
the X-ray luminosity threshold in the definition of an X-ray fossil, 
scatter in Lx will alter the X-ray fossil number density, since the 
number density of haloes is a steep function of mass. 

• By selecting optical fossils from groups randomly generated 
from a Schechter luminosity function, we demonstrate that for 
small numbers of galaxies per group, a significant fraction of opti- 
cal fossil groups are expected to be purely statistical, requiring no 
special physical mechanism to generate the 2 magnitude luminos- 
ity gap. For groups with more than 40 members, this effect largely 
disappears, with very few fossil groups expected at random. 

• The probability of finding optical fossils with mass Ad, i.e, 
P/(Af; optical) is a decreasing function of: (a) group dark matter 
halo mass and, (b) the fraction of the virial radius within which the 
first and second brightest galaxies are being found. Conversely, as 
dark matter halo mass become small, the probability Pf (M\ X-ray) 
for X-ray fossils decreases. 

• Both high-mass and low-mass X-ray fossil groups are found 
to have assembled ~ 90% of their final masses by a redshift of 
2 = 0.24. The corresponding mass fraction is about 70 — 80% 
for two different sets of high-mass control samples, and ~ 85% 
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Figure 7. Evolution of a typical of massive X-ray fossil group (right) in 
comparison to a typical massive normal group (left) from redshift 2=1.0 
to 0. The dimension of each panel is 10 X 10 Mpc. The poi nts represent 
indivi dual gas particles from the Millennium Gas simulation jPearce et alj 
I2OO7I) . 



for low-mass control samples, where groups fulfill the same X- 
ray luminosity criterion (> 3.34 x 10^* /i^^M©) but have the 
optical luminosity gaps corresponding 0.1 ^ Ami2 ^ 0.3 and 
0.8 ^ Ami2 1.0 magnitudes. 

This study shows that fossils indeed are formed early, with 
more than ~80% of their mass accumulated as early as 4 Gyr 
ago. They are also relatively isolated compared to non-fossils. The 
strongest X-ray fossil candidates are those with the highest X-ray 
luminosity as these systems are not expected to have a large lumi- 
nosity gap between their first and second ranked galaxies entirely 
by chance. As always, systems with more than a handful of galaxies 
are to be preferred. 

In principle, comparison of the observed space density of fos- 
sils as a function of X-ray luminosity with that of fossils from sim- 
ulations can provide valuable constraints on the treatment of phys- 
ical processes included in the simulations. The tentative evidence 
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Figure 8. Comparison between the mass growth of high and low-mass X- 
ray fossils plotted in Figs.|6j) and |6j;. 



for a discrepancy, whereby the observed space density of low X- 
ray luminosity fossil groups may exceed that predicted from the 
simulations, will be worth revisiting in the future, when better ob- 
servational estimates are available. 

It is interesting that while the amount of recent mass assem- 
bly in control groups increases with halo mass, as is expected in 
the hierarchical growth paradigm, there is almost no difference be- 
tween the mass assembly of high-mass and low-mass X-ray fossils 
after redshift z = 0.6 (see Fig.[8](. It seems that both low-mass and 
high-mass fossil systems are undisturbed at low redshift. 

Since we expect faster orbital decay and more efficient galaxy 
merging in lower mass systems, due to the lower velocity disper- 
sion of individual galaxies within the group, we would expect to 
find a higher incidence rate of fossils amongst poor groups. Fig[3] 
shows that this is indeed the case, but the effect is not very strong, 
once the influence of statistical fossils is removed, and in X-ray 
fossils any rising trend at low richness is overwhelmed by the fact 
that many of the optical fossils fail to exceed the X-ray luminosity 
threshold. A word of caution about the treatment of orbital decay 
is in order here. As discussed in Sec. 12.21 the orbital evolution of 
subhaloes in the Millennium simulation is well treated until these 
subhaloes are reduced by stripping to 20 dark matter particles, but 
thereafter is calculated semi-analytically, using an approximate for- 
mula. In practice, a significant fraction of the second ranked galax- 
ies in the fossil systems we have extracted from the simulation have 
been stripped below this 20 particle limit. For example, in fossils 
with masses of only 10^^ Mq, approximately 35% of second 
ranked galaxies have been stripped below the limit, though this 
fraction drops to ~I0% for systems with mass > 10^* h^^ Mq. 
For such galaxies the timescale for their subsequent decay and 
merger with the central galaxy is not very reliable. However, this 
is not a major issue for massive haloes, and it is interesting and sur- 
prising that the incidence of fossils in rich systems is fairly flat at 3- 



4%. Observational studies should, in due course, show whether this 
is reflected in the real Universe. On e example of a fairly rich f ossil 
cluster has already been reported bv Khos roshahi et"alT ( l200d) . 

The magnitude gap distribution of haloes at different X-ray lu- 
minosities and the mass evolution of fossil groups discussed above 
both support the idea that X-ray fossil groups are not a distinct 
class of objects but rather that they are extreme examples of groups 
which collapse early and experience little recent growth, so that 
their galaxies have time to undergo orbital decay and merging. The 
X-ray and optical scaling properties of such extreme groups can be 
expected to differ from those of groups with more typical evolu- 
tionary histories, and such differences have already been observed 
jKhosroshahi. Poimian & Jonesl l2007). A comparison of such ob- 
served differences with the properties seen in the Millennium sim- 
ulation groups is underway, and should provide a valuable check 
on the adequacy with which feedback processes and other baryon 
physics is handled in the simulations. 
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